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EXTREME  VALUES  OF  BIRTH  AND  DEATH  PROCESSES  AND  QUEUES 

by 

Richard  F.  Serfozo 
Georgia  Institute  of  Technology 

Abstract 

We  study  the  asymptotic  behavior  of  maximum  values  of  birth  and 
death  processes  over  large  time  intervals.  In  most  cases,  the 
distributions  of  these  maxima,  under  standard  linear  normalizations, 
either  do  not  converge  or  they  converge  to  a  degenerate  distribution. 
However,  by  allowing  the  birth  and  death  rates  to  vary  in  a  certain 
manner  as  the  time  interval  increases,  we  show  that  the  maxima  do  indeed 
have  three  possible  limit  distributions.  Two  of  these  are  classical 
extreme  value  distributions  and  the  third  one  is  a  new  distribution. 

This  third  distribution  is  the  best  one  for  practical  applications.  Our 
results  are  for  transient  as  well  as  recurrent  birth  and  death  processes 
and  related  queues.  For  transient  processes,  the  focus  is  on  the  maxima 
conditioned  that  they  are  finite. 


Keywords:  Extreme  values,  birth  and  death  processes,  M/M/s  queues,  limit 

|  F 

theorems . 

;  i 

} 

This  research  was  sponsored  in  part  by  grant  AFOSR  84-0367 J 


?0R Cl  OFFICE  OF  SCIENTIFIC  RESEARCH  (AFSC) 
CE  OF  TRANSMITTAL  TO  DII C  and  lg 

:ribution  is  unlimited. 

;  Information  Division 


T;:. 


ApprovJd^orptbUorslease 

distribution  unli»it-*d* 


1.  Introduction 


When  modeling  the  dynamics  of  a  parameter  of  a  system  by  a 
stochastic  process,  the  questions  one  addresses  depend  on  the  nature  of 
the  parameter.  In  some  instances,  the  extreme  values  of  the  parameter 
rather  than  its  usual  values  may  be  of  paramount  interest.  In  a 
manufacturing  plant,  for  example,  a  typical  parameter  is  the  queue  length 
of  parts  waiting  to  be  processed  at  a  work  station.  Small  to  moderate 
values  of  the  queue  may  indicate  that  the  system  is  operating 
successfully  and  the  queue  fluctuations  are  unimportant.  On  the  other 
hand,  large  queues  may  call  for  extraordinary  measures  such  as  allocation 
of  auxiliary  storage  space,  employee  overtime,  or  rescheduling  of 
production.  A  natural  question  is:  What  is  the  probability  that  the 
queue  will  exceed  a  specific  critical  value  in  a  certain  time  period? 
Extreme  value  questions  like  this  are  the  topic  of  this  paper.  More 
specifically,  our  focus  is  on  characterizing  the  asymptotic  behavior  of 
the  maxima  of  birth  and  death  processes  and  related  queues. 

The  gist  of  our  study  is  illustrated  by  the  following  results  for 
the  M/M/s  queue.  Consider  such  a  queueing  process  in  which  customers 
arrive  to  s  servers  according  to  a  Poisson  process  with  rate  X,  and  the 
independent,  exponentially  distributed  service  times  have  mean  p  *.  Let 
Mn  denote  the  maximum  queue  length  in  the  time  interval  10, T^],  where 
is  the  nth  time  the  system  becomes  empty.  Our  interest  is  in  finding 
norming  constants  a^,  bfl  >  0  and  a  non-degenerate  distribution  G  such 
that 

(1.1)  lim  P((M  -  a  )/b  <  x)  -  G(x), 

n—  n  n  n 

for  each  continuity  point  of  G.  When  such  a  ,  b  ,  G  exist,  we  say  that 


i 


! 


M_  has  the  limit  distribution  G.  Otherwise,  we  say  that  M  does  not  have 
n  - - n - 

a  limit  distribution.  As  usual,  we  consider  only  linear  normalizations. 

We  can  write  M  *  max{Y,,...,Y  r  ,  where  Y,  is  the  maximim  of  the 
n  I  n  k 

queue  in  the  interval  [T^pT^l*  Since  the  queuing  process  is  Markovian, 
then  Y^Y^,..*  are  independent  identically  distributed  random  variables. 
From  the  classical  extreme  value  theory  for  independent  identically 
distributed  variables  (see  for  instance  Galambos  (1978)  or  Leadbetter  et 
al.  (1983)),  we  know  that  the  possible  limit  distributions  for  are 

•  v 

only  exp(-x  ),  x  >0,  or  exp(-e  ),  xeR.  One  consequence  is  as  follows; 
this  is  a  special  case  of  Theorem  2.5. 

THEOREM  1.1.  If  the  queueing  process  is  null  recurrent  (A=  sp),  then 

_  -1 

lim  P(M  /(n(  A/y)Ss!)  <  x)  ■  e~x  ,  x  >  0. 
n+« 

Cohen  (1969)  proved  an  analogue  of  this  for  the  M/G/l  and  G/M/l  queues; 
related  studies  on  extreme  values  of  queue  lengths  and  waiting  times  are 
Heyde  (1971)  and  Iglehart  (1972). 

Another  classical  result  for  discrete  random  variables  is  that  the 
convergence  (1.1)  can  take  place  only  if 

lim  P(Y.  »  m)/P(Y1  >  m-1)  *  0. 

Consequently,  we  have  the  following  anomaly. 

THEOREM  1.2.  If  the  queueing  process  is  positive  recurrent  (A<  sp), 

then  M  does  not  have  a  limit  distribution, 
n 

This  non-convergence  theorem  is  surprising,  especially  in  light  of 
Theorem  1.1,  since  positive  recurrent  processes  generally  have  nicer 
properties  than  null  recurrent  ones.  Cohen  (1969)  and  Anderson  (1970) 


give  insights  on  typical  liminf's  and  limsup' s  for  the  distribution  of 


(M  -a  )/b  . 
n  n  n 


Our  point  of  departure  is  to  establish  the  convergence  (1.1)  for 
queues  and  birth  and  death  processes  in  spite  of  the  non-convergence 


described  in  Theorem  1.2  and  its  generalization,  Theorem  2.3.  Our 


approach  is  to  allow  the  birth  and  death  parameters  (X,p,s  for  the  M/M/s 

queue)  to  vary  with  n  when  considering  the  convergence  of  M  .  Such 

n 

parameter  variations  in  limit  theorems  are  not  uncommon.  A  classic 


example  is  that  if  is  a  binomial  random  variable  with  parameters  n,p 

and  if  p  =  p  varies  with  n  such  that  np  ♦  1  >  0  as  n  ♦  «,  then  the 
n  n 

d istribut ion  of  converges  weakly  to  a  Poisson  distribution  with  mean 

X. 


Here  is  an  example  of  our  major  results.  Suppose  Mn  is  the  maximum, 

as  above,  of  an  M/M/s  queue,  where  the  arrival  rate  X  -  X(n)  and  service 

rate  p  =>  p(n)  depend  on  n.  Let  pn  *  X(n)/(sy(n)) . 

THEOREM  1.3.  Suppose  that  p  <1  for  each  n  and  that  p  ♦  1.  The 

n  n 

possible  limit  distributions  for  M  are  G_(x)  =  exp(-x-1),  x  >  0;  G  (x) 

n  u  ® 

=  exp(-e  ),  xeR;  and 

(1.2)  GQ(x)  “  exp(-c/(eX-l))  ,  x  >  0,  for  0  <  c  <  ». 

The  has  the  limit  distribution  Gc ,  where  0  <  c  <  <■,  if  and  only  if 

n(l-p  )/s!  ♦  c.  Appropriate  nonning  constants  a  ,  b  are  as  in  Theorem 
n  n  n 

3.1. 

Note  that  GQ  and  Gw  are  classical  extreme  value  distributions.  The 
third  distribution  (1.2)  has  not  appeared  in  the  literature  before. 
Numerical  work  has  shown  that  this  new  distribution  is  the  best  one  for 
practical  approximations.  Namely,  for  the  standard  M/M/s  queue  with 
traffic  Intensity  p»  X/sp, 


4 


a 


s 

ft 


P<-Mnlogp  <  x)  s  exp(cn/(e  -i)),  x  >0, 

where  c  *  n(l-p).  This  approximation  is  good  for  n  >  15  and  for  all 
n 

P  <  1. 

This  completes  our  introduction.  Here  is  what  lies  ahead:  Section 
2  consists  of  preliminaries,  including  classical  convergence  and 
non-convergence  theorems  for  extreme  values  of  birth  and  death  processes; 
Section  3  contains  our  main  results  for  recurrent  birth  and  death 
processes;  Section  4  contains  analogous  results  for  transient  processes; 
and  Section  5  gives  applications  to  M/M/s  and  related  queues. 


2.  Preliminaries 

We  shall  consider  a  continuous-time  birth  and  death  process  on  the 

nonnegative  integers  with  birth  rates  X^.X^.X^,...  a°d  death  rates 

Mq  ”  0,  Uj,^****  This  is  a  Markov  process  that  evolves  as  follows: 

Upon  entering  state  k,  the  process  remains  there  for  an  exponentially 

distributed  time  with  mean  (X^+p^)  * ,  and  then  it  moves  to  state  k+1  or 

k-1  according  to  the  respective  probabilities  X^/(  X^+p^)  and  ^/(X^M.)* 

We  assume  for  now  that  the  process  is  positive  recurrent  or  null 

recurrent  -  transient  processes  are  discussed  in  Section  4.  We  also 

assume,  for  convenience,  that  the  process  at  time  zero  begins  in  state 

zero.  Let  M  denote  the  maximum  value  of  the  process  in  the  time 
n 

interval  [0,T  ],  where  T  is  the  time  of  the  nth  visit  of  the  process  to 
n  n 

state  zero.  The  Mn  and  Tn  are  finite  valued  since  the  process  is 

recurrent.  We  shall  study  the  asymptotic  behavior  of  the  distribution  of 

M  as  n  ♦  0. 
n 

We  can  write  Mn  *■  max  {Y^  , . . .  ,Yn },  where  Y^  is  the  maximum  of  the 
birth  and  death  process  in  the  time  interval  [T,  ,T,],  here  Tn  -  0. 

K"  1  1C  U 

Since  the  process  is  Markovian,  the  random  variables  Y  ,Y  ,...  are 

1  2 


.  < . 

■r.  '-.V. 


5 


m 


|«S 

i*S 


i 


p 


■Si 


independent  and  identically  distributed.  Consequently, 

P(Mn  <  x)  =  P(Yj  <  x,...,Yn  <  x)  =  F(x)", 
where  F  is  the  distribution  of  the  Y^'s.  To  obtain  an  expression  for  F, 
note  that  the  successive  states  of  the  birth  and  death  process  form  a 
simple  discrete-time  random  walk  that  moves  from  state  k  to  state  k+1  or 
k-1  according  to  the  respective  probabilities  X^/(X^+p^)  and  n  /(\^+p^). 
Then  clearly  F(x)  is  the  probability  that,  starting  from  state  1,  the 
random  walk  reaches  state  0  before  it  exceeds  x.  Thus,  from  Section  1.12 
of  Chung  (1967),  we  know  that 


(2.1) 


F(x)  =  1  -  (  Z  r  )‘ 
k-0 


where  r^  =  1  and  r^  ■  (  . .  UJt)/(  X^  •  •  •  X^)  ,  k  >  1.  Furthermore,  the  birth 

00 

and  death  process  is  recurrent  when  ^jc=grj(  *  00  and  is  transient  when 


Z.  <  °°« 
k=0  k 


Here  are  some  asymptotic  properties  of  the  ratio  X^/p^,  depending  on 
whether  the  birth  and  death  process  is  transient  or  recurrent;  we  let 


(2.2) 


p  ■  liminf  X,/  and  p  *  limsup  X^/  p^. 
k-*»  k+» 


LEMMA  2.1.  If  <  *,  then  X^/  >  1  for  an  infinite  number  of  k' s , 

and  p  >  1.  If  *  •»  then  p  <  1. 

Proof.  Suppose  r^  <  «.  For  any  N  >  1 ,  we  have 

00  00 

rN-l  2  ^  ‘‘N*  *  *  *N*  *  *  *  2  rk  <  "* 

N  1  k=N  k=N 

Thus  X^/ p^  must  be  >  1  for  an  infinite  number  of  k's;  otherwise,  the 
first  sum  would  be  infinite.  To  prove  p  >  1,  fix  e  >  0  and  let  N  be  such 


■ 


that  X^/ >  p  -  e,  k  >  N.  Then 

os  N-l  • 

•  >Er  »  Er  +r  E(p  -  e)_m. 
k=0  k=0  m=0 

Consequently,  ( p- e)  ^  <  1  for  each  e  >  0,  and  hence  p  >  1.  The  second 
assertion  Is  proved  similarly. 

Keep  in  mind  that,  for  now,  we  are  considering  only  recurrent  birth 
and  death  processes.  Our  interest  is  in  the  weak  convergence  of  the 
distribution 

(2.3)  P( (M  -  a  )/b  <  x)  -  F(an  +  b  x)n  -  [1  -  (1  -  FU  +  b  x))]", 

nnn  nn  nn 

where  a  and  b  >  0  are  constants.  It  is  well  known  that,  for  any  y  eR 
n  n  n 

and  -<*>  <  y  <  •,  the  convergence  (1+y  )n  ♦  eY  is  equivalent  to  ny  ♦  y. 

n  n 

This  property  applied  to  (2.3)  translates  into  the  following  known  result 
(cf  Corollary  1.3.1  of  Galambos  (1978)  or  Theorem  1.5.1  of  Leadbetter  et 
al.  (1983)).  Here  0  <  g(x)  <  ». 

CONVERGENCE  CRITERION  2.2.  P((M  -a)/b  <  x)  e"B(x>  as  n  ♦  »  if  and 

n  n  n 

only  if  n(l  -  F(afi  +  b^x))  ♦  g(x)  as  n  ♦  •».  Note  that  this  criterion  is 
also  true  when  F  varies  with  n. 

We  begin  by  characterizing  when  the  maximum  does  not  have  a  limit 
distribution  and  when  it  might  have  one.  Here  we  use  the  notation  (2.2) 


^  *=k+l  *  1 

THEOREM  2.3.  (i)  The  maximum  M  does  not  have  a  limit  distribution  if 

n 


and  only  if 


(2.4) 


m-1 

liminf  E  a  <  «. 
.  ~  ran 


CP  •m 

(ii)  Inequality  (2. A)  holds  when  E^q  lirasup  <  «  or  when  p  <  1, 


(iii)  Inequality  (2.4)  does  not  hold  when  liminf 

m  -**> 


®  or  when 


p  >  1.  The  condition  p  >  1  is  equivalent  to  lim  X  /a  =  1. 

k+«  u  k 


(iv)  If  =  lim  exists  for  each  k,  then  inequality  (2.4)  holds  if 
m+® 

00 


and  only  if  Ek=0  \ 
Proof.  ( i)  We  know 
that  there  exist  x 

n 

and  only  if 


<  ®. 

,  from  Theorem  1.7.13  of  Leadbetter  et  aU  (1983), 
and  0  <  t  <  ®  such  that  n(l  -  F(x  ))  +  t  as  n  +  •  if 


(F(m)  -  F(m  -  1 ) ) /( 1  -  F(m  -  1))  +  0  asm  •*■<*>. 
The  latter  condition,  in  light  of  (2.1),  is 


l  m~*  l 

VCo  rk  '  : £„  “kj"  *° 

k=0 

Then,  by  Criterion  2.2,  the  M  does  not  have  a  limit  distribution  if  and 

only  if  the  last  convergence  does  not  hold,  which  is  equivalent  to  (2.4). 
00 

(ii)  If  the  sum  £j_  limsup  is  finite,  then  by  Fatou' s  lemma  this 

n-H» 


finite  sum  is  an  upper  bound  for  limsup 

m+® 


jf"1 

k-0 


akm 


and  hence  (2.4)  holds. 


Now  suppose  p  <  1. 
number  N  such  that 

(2.5) 


Then  for  any  c  in  the  interval  (p,l), 


X^/ <  c  for  k  >  N.  It  follows  that, 

N  m 

E 


-1 


£o.=r  E  r.  +  j,  a. 

,  _ran  m  .  n  k  .  tan 

k=0  k=0  k=N+l 

,  „  N  m-N 

<  ^  cm'N  E  r  ♦  £  c* 

N  k=0  t =0 


+  l/(  1  -  c)  as  ra  ♦  ®. 


there  is  a 
for  m  >  N, 


Thus  (2.4)  holds. 

(iii)  By  Fatou' s  lemma,  the  sura  E^q  liminf  is  a  lower  bound  for  the 

m+® 

left  side  of  (2.4).  Thus,  if  this  sum  is  infinite,  then  (2.4)  does  not 


hold.  Now  suppose  p  >  1.  Then  for  any  c  <  1,  there  is  an  N  such  that 
X^/ >  c  for  k  >  N.  The  display  (2.5)  clearly  holds  with  the  inequality 


reversed  and  so 


liminf  £ 
k=0 


>  l/(l-c). 


Letting  c  1  implies  that  (2.4)  does  not  hold.  To  prove  that  p  >  1  if 

and  only  if  lim  X,  /  =  1 ,  we  need  only  show  that  p  >  1  implies 

k-*-» 

p  *  p  m  1.  But  this  is  true  because,  by  Lemma  2.1,  we  know  that 
p  <  p  <  1. 

(iv)  This  part  is  a  consequence  of  parts  (ii)  and  (iii). 

The  preceding  theorem  yields  the  negative  conclusion  that  the 
existence  of  a  limit  distribution  for  is  the  exception  rather  than  the 
rule:  There  is  no  limit  distribution  for  a  typical  positive  recurrent 

process  with  lirasup  X^/ <  1,  but  one  might  exist  for  an  atypical 
k-M»  S 

process  with  lim  X^/ =  1  (such  as  a  null  recurrent  process).  For  those 
k-H* 

instances  when  there  might  be  a  limit  distribution  for  M^,  we  have  the 

following  properties  from  the  classical  extreme  value  theory. 

PROPERTIES  OF  M  2.4.  (a)  The  possible  limit  distributions  for  M  are 
n  n 

—  y 

only  exp  (-x  ),  x  >  0,  or  exp  (-e  ),  xeR. 

(b)  The  first  of  these  distributions  is  the  limit  if  and  only  if 

t  tx 

(2.6)  (1  -  F(  tx) )  /( 1  -  F(  t) )  =  £  r  /  £  r.  -*•  xY  as  x-*-®. 

k=0  k=0 

Appropriate  norming  constants  are  a^  =»  0  and  bn  =  min{m:  1^  r^  >  n}. 


(c)  The  second  distribution  in  (a)  is  the  limit  distribution  for  M  if 

n 

and  only  if  there  is  a  positive  function  g(t)  such  that 


(2.7) 


t  t+xg(t) 

(l  -  F(t  +  xg(t)))/(l  -  F(  t) )  =  £  r  /  £ 

k=0  k=0 

-x 

♦  c  as  t  +  «. 


In  this  case,  one  can  choose 


oo  k 


g<t> "  .V*  „V  £nt»)  • 

k=0  k=t  m=0 


and  a  =  min{m:  £T  „r.  >  n)  and  b  =  g(a  ). 

n  k=0  k  n  n 

A  special  case  of  ( b)  is  as  follows.  This  applies,  for  instance,  to 
a  null  recurrent  process  with  =  u^,  k  >  s,  for  some  s  (such  as  the 
M/M/s  queue  in  Theorem  1.1). 


THEOREM 


2.5.  If  £=0  |1  -  yuj  <  »,  then 


(2.8) 


lim  P(M  /nb  <  x)  =  e  ,  x  >  0 
n 

n-*» 


where  b  3  II  X,/  a  . 
k=l 


Proof.  A  basic  property  of  infinite  products  of  real  numbers  is  that 


II  (l-a^)  exists  when  <  *•  In  light  of  this,  the  hypothesis 

k=l 

implies  the  existence  of  the  limit  b.  Now 


lira  r  =  lim  II  p  /  A  =  b  . 
k+«  K  k-~°  fc=l  *  X 


Consequently,  n  *  £^_q  rk  +  b  ^  *  anc*  so 


-l  nbx  _i  _i 

n(l  -  F(nbx)  ]  =  [bx(nbx)  £  r  ]  ♦  x"  as  n  > 


This  convergence  and  Criterion  2.2  yield  (2.8). 


3.  Main  Results 


We  saw  above  that  the  maximum  M  does  not  have  a  limit  distribution 

tl 


for  a  wide  class  of  birth  and  death  processes,  including  those  in  which 


limsup  X^/ p^  <1.  However,  might  have  a  limit  distribution  when 
k-H» 

lim  \  =1.  These  negative  and  slightly  positive  findings  prompted 
k-*» 

us  to  explore  the  convergence  of  when  the  parameters  X^  =  X  and 
p  =  p  vary  with  n  such  that  X  /  p  is  nearly  1  for  large  n  and  k. 

K  UK  Iliv  uK 

This  is  the  basis  of  our  following  results. 

Consider  a  sequence  of  recurrent  birth  and  death  processes  indexed 

by  n  =  1,2,...,  where  the  nth  process  has  the  respective  birth  and  death 

rates  X^  and  p^  when  in  state  k  =  0,1,...  For  the  nth  process,  let 

denote  its  maximum  up  to  the  time  of  its  nth  return  to  state  zero.  This 

M  has  the  same  meaning  as  the  one  in  Section  2,  but  here  its  defining 
n 

parameters  X  ,  ,  p  ,  vary  with  n  as  well  as  with  k.  That  is,  M  is  the 
nk  nk  n 

maximum  of  n  independent  random  variables  with  the  common  distribution 

F  (x)  -  1  -  (  Z  r  )_1 
n  k=0  * 

where  r  „  =  0  and  r,  =  (p  ,  . . .  p  .  )/(  X  , . . .  X  ,  )  ,  k  >  1. 
nO  nk  nl  nk  nl  nk 

We  shall  assume  that,  for  each  n,  there  is  a  positive  number  p  <  1 

n 

and  a  positive  integer  s^  such  that 


(3.1) 


X  ,  /  p  ,  =  p  for  k  >  s  , 
nk  nk  n  n 


(3.2) 


s 

lim  n_1  Zn  r  .  =  0. 

i  n  nk 

n  -**>  k=0 


The  parameter  s  may  be  bounded,  unbounded,  or  independent  of  n,  such  as 
n 

oo  oo 

s  =1.  Assumption  (3.1),  with  p  <  1,  implies  that  L,  ^r  ,  >  E,  p 

n  ’  n  k=0  nk  k=s  *n 


=*  ",  which  ensures  that  the  nth  birth  and  death  process  is  recurrent. 
Assumption  (3.1)  is  satisfied  automatically  for  M/M/s  queues:  the  s  is 


the  number  of  servers.  Assumption  (3.2)  holds  when  s  and  r  are 

n  nk 


bounded  or  when  r  ,  <  B  for  k  <  s  and  s  B  /n  ■*•  0. 

nk  n  n  n  n 


We  shall  show  that  when  -►  1 ,  the  possible  limit  distributions  for 


(M  -  a  )/b  are  as  follows.  Here  c  =  lim  c  and 
n  n  n  n 

n-*-" 


s 


(3.3) 

Case 


c  =  n(l  -  p Jr"1  =  n(l  -  p  ),  n”  X  .  /  y  .  . 
n  n  ns  n  k=l  nk  nk 

n 


Limit  Distribution 


Norming  Constants 


c  =  0  Gf,(x)  =  exp(-x  *),  xX>,  a  =  s  -1  b  “nil11!,/ 

U  n  n  n  ,  .  nk 


k=l 


nk 


0  <  c  <  "  G  (x)  =  exp(-c/(e  -1))  a  =s  -  1  b  =  -l/logp 
c  n  n  n  n 


C=  oo 


-J( 

G  (x)  =  exp(-e  ),  xeR,  a  =  s  -  1  -logc  /logp 
00  nn  n  n 


b„  '  -/l08V 


An  easy  check  shows  that  the  distributions  G  (0  <  c  <  ®)  are  of  distinct 

c 

type:  Gc  and  Gc,  are  of  the  same  type  (G£(x)  =  Gc,(a  +  bx)  for  some  a,b) 
if  and  only  if  c  ■  c'. 

The  following  result  gives  sufficient  conditions  for  the  existence 

of  limit  distributions  for  M  . 

n 

THEOREM  3.1  Suppose  (3.1)  and  (3.2)  hold  and  p  +1.  If  c  ♦  c  as 

n  n 

n  where  0  <  c  <  then 

(3.4)  lim  P  f(M  -  a  )/b  <  x }  =  G  (x),  xeR, 

b  n  n  n  '  c 

n+® 

where  the  Gc>  a^,  b^  corresponding  to  the  limit  c  are  displayed  above. 

Proof.  By  the  Convergence  Criterion  2.2,  it  suffices  to  show  that  c  ♦  c 

n 

implies 


lim  n(l-F  (a  +  b  x))  =  x 
n  n  n 

n  +« 

*  c/(eX-l) 


when  c  =  0 
when  0  <  c  <  " 


(3.5) 


=  e  when  c  *  <*>. 

To  this  end,  let  m  (x)  denote  the  integer  part  of  a  +  b  x.  Usin; 

n  n  n 


(3.2),  we  have 


(3.6) 


_i  mr/x) 

n( 1  -  F  (a  +  b  x))  -  (n  Z  r  ) 
n  n  n  ,  _  nk 

k=0 

,  s  -1  ,  m  (x)  v  , 

r  -1  _n  ,  -1  n_  -( k-s  )  v-1 

=  In  r  r,+nr  £  p  n  I 

.  nk  ns  ,  Kn  ' 

k=0  n  k=s 

n 

■  (o(l)  +  zn(x)) 


where 


s  -m  (x) 

z(x)  =  (p  n  n  -p)/c  when  p  <  1 
n  n  n  n  n 

=  (n  (x)  -  s  )/(nr  *)  when  p  =1. 

n  n  ns  n 

n 

Then  to  establish  (3.5),  it  suffices  to  show  that  zn(x)  *  converges  to 


the  values  on  the  right  side  of  (3.5).  Three  cases  present  themselves 

Case  1 :  c  *  c  ■  0.  Here 
n 

m  (x)  -  a  +  b  x  +  0(1)  -  s  -  1  +  xnr  1  +  0(1). 

n  n  n  n  ns 

n 

First  consider  the  special  situation  in  which  *  1  for  each  n.  Usin 

n  *r  +  0,  we  have 
ns 

n 

z(x)  1  =  nr  1  /(m(x)-s) 
n  ns  n  n 

n 

*  (x  +  o(l))  *  +  x  *  as  n-H». 

Next  consider  the  general  situation  in  which  p^  -*■  1.  Because  of  the 
preceding,  we  may  assume  that  p^  <  1  for  each  n. 

Using 

log  pn  =  -  (1  -  Pn)  +  o(l  -  pn)  and  eU  =*  1  +  u  +  o(u) 


as  u  >  0,  it  follows  that 


s  -m  (x) 

pn°  n  '  pn  =  expf(sn  "  nn(x))logpJ  ~  P„ 

=  exp[xc  +  0(1 —  p  )]  -  p 
n  n  n 

=  1  +  xc  +  0(l-p  )  +  o(c  )  -  p 
n  n  n  n 

=  xc  +  0(l-p  )  +  o(c  ). 
n  n  n 

Substituting  this  in  the  expression  for  z  (x),  and  using  (1-p  )/c 

n  n  n 

■  n  *r  ♦  0,  we  have  the  desired  convergence 
n 

z  (x)  *  =  [x  +  0 ( 1—  p  )/c  +  o(l)]  *  -*■  x  *  as  n  +  ®. 

n  n  n 

Case  2 :  c  -*■  c  and  0  <  c  <  “.  Here 
n 

m  (x)  =  s  -  1  -  x/logp  +  0(1). 
n  n  n 


Using  logp^  53  o(l),  we  have 


zn(x)  =  cn/(exp[(sn  -  mn(x))logpn]  ~Pn) 

=  c  /(exptx  +  o(l)]  -  p  ) 
n  n 

♦  c/(e  -  1)  as  n 


Case  3:  c  *  c  «  ®.  Here 
n 


m  (x)  *  s  -  1  -  x/logp  +  logc  /logp  +  0(1). 
n  n  n  n  n 


Using  logp  =  o(l),  we  have 
n 


s-m(x) 

P  n  "  =  exp{  ( s  -  m  (x)  )logp  ]  -  c  ex  \ 

n  n  n  n  n 


,  x-1  ,  x+o(  1 )  ,  .-1  -x 

z  (x)  =  (e  -  p  / c  )  ♦  e  as  n  ♦  ®  , 

n  n  n 

This  completes  the  proof. 

Theorem  3.1  says  that  the  convergence  of  c^  is  sufficient  for  to 

have  a  limit  distribution  when  p  +  1.  The  next  result  says  that  the 

convergence  of  c  is  necessary  as  well,  and  that  M  has  no  limit 
n  n 


few 


distribution  other  than  those  above. 

THEOREM  3.2.  Suppose  (3.1)  and  (3.2)  hold  and  ♦  1.  The  possible 

limit  distributions  for  M  are  G  ,  0  <  c  <  «;  and  M  has  the  limit 

n  c  n 

distribution  G  if  and  only  if  c  ♦  c  as  n  -►  °». 
c  n 

Proof.  Suppose  has  the  limit  distribution  H.  Since  [0,°»]  is  a  closed 

set  in  the  extended  real  line,  there  are  positive  integers  +  <*>  and  c 

in  [0,®]  such  that  c  +  c  as  k  >  ®.  Then  by  Theorem  3.1,  we  know  that 

nk 

M  has  the  limit  distribution  G  .  Moreover,  M  also  has  the  limit 
nk  c  nk 

distribution  H.  From  Khintchine's  theorem  on  convergence  to  types  of 

distributions  (see  for  instance  Theorem  1.2.3  of  Leadbetter  et  al. 

(1983)),  it  follows  that  H  and  G  are  of  the  same  type.  This  proves  that 

any  limit  distribution  of  M  must  be  one  of  the  distributions  G  , 

n  c 

0  <  c  <  “>. 

We  now  prove  that  c^  ♦  c  is  necessary  and  sufficient  for  to  have 

a  limit  distribution.  The  sufficiency  follows  from  Theorem  3.1.  To 

prove  the  necessity,  suppose  that  has  the  limit  distribution  Gc«  Let 

c  be  any  convergent  subsequence  of  c  and  let  c'  *  lim  c  .  Arguing  as 
\  n  k+~  \ 

in  the  last  paragraph,  it  follows  that  Gc,  as  well  as  Gc  is  a  limit 

distribution  of  M  and  that  G  ,  and  G  are  of  the  same  type.  Conse- 
n  c’  c 

quently,  c'  ■  c.  Thus,  we  have  shown  that  any  convergent  subsequence  of 


I 


m 


I 


m 


'f'j 
V -J 

J-., 


c^  must  converge  to  c  and  hence  c^  ♦  c. 


APPROXIMATION  OF  P(M  <  x)  FOR  PRACTICAL  APPLICATIONS  3.3.  Consider  the 

n 


maximum  of  a  birth  and  death  process  with  rates  (without  the 

artificial  dependence  on  n)  that  satisfy  ■  p,  k  >  s,  for  some  0  <  p 

<  1  and  s  >  l.  Theorems  3.1  and  3.2  yield  the  approximation 


(3.7)  P((M  -  s  +  l)(-logp)  <  x)  2  exp(c  /(e  -  1)),  x  >  0, 

n  n 


.**  V 


v> 


5 


where  c  =  n(l-p)  II  and  n  is  large. 

k»l 

There  are  analogous  approximations  for  P(M  <  x)  by  G_  when  c  is 

n  On 

small  and  by  Ga  when  c^  is  large.  However,  these  are  not  as  good  as 

(3.7),  which  is  superior  for  any  c^.  This  is  because  and  Gw  are  only 

theoretical  limits  for  the  two  "unobtainable"  values  of  c^  in  [ 0 , °°} ;  they 

are  not  functions  of  the  actual  c^  as  the  right  side  of  (3.7)  is.  For 

the  case  when  the  process  is  null  recurrent  ( p=  1),  Theorem  2.4  yields  the 

s 

classical  approximation  P(Mn/nb  <  x)  =  Gq(x),  where  b  *  n  X^/m^*  We 

k*l 

were  pleasantly  surprised  that  the  approximation  (3.7)  is  accurate  even 
when  p  is  not  near  one.  This  is  apparently  because  p  appears  on  the 

right  as  well  as  left  side  of  (3.7).  For  the  M/M/1  queue,  we  found  from 

numerical  computations  that  the  difference  between  the  two  sides  of  (3.7) 
is  below  0.018  when  n  »  15  and  below  0.01  when  n  **  20,  for  any  p  in 
(0,1). 

PROPERTIES  OF  THE  LIMIT  DISTRIBUTIONS  3.4.  Let  X  denote  a  random 

c 

variable  with  distribution  Gc>  0  <  c  <  and  let  Y  denote  an 

exponentially  distributed  random  variable  with  unit  mean.  Standard 

chang e-of-variable  computations  show  that 

XQ  ^  i/Y,  -logY,  Xc  2  log(Y  +  c)  -  logY,  0  <  c  <  «, 

where  these  are  equalities  in  distribution.  Solving  the  equations  for  Y 
and  using  obvious  substitutions,  we  also  have 


XQ  =*  =  exp(-Xco)  2  c/(exp(X£)  -  1) 

X  £  logX-  log((exp(X  )  -  l)/c) 

00  U  c 

X  logo  +  cXfl)  log(l  +  cexp(-X<J). 


The  G  can  be  viewed  as  the  limit  of  G  as  c  +  «•  in  that 

•  c 


m 


G  (x  +  logc)  ♦  G  (x)  as  c  ♦  •  (i.e.  X  -  logc  ►  X  )•  It  Is  known  that 
c  09  c  * 

oo 

EX m  ■  E(-logY)  ■  -  /  e  y  logy  dy  »  y, 

0 

ra 

where  y  =>  .5772...  Is  the  Euler-Mascheroni  constant  y  -  lim  (  E  1/n  - 

m  n=  1 


logra)  ,  and 


ECe0*-)  =*  E(  Y~ ^  -  r(l-a),  for  a<  1, 


where  T(  a)  *  /  u0*"^  u  du.  In  comparison,  for  0  <  c  <  »,  we  have 

0 

EX  »  E( log( c  +  Y)  -  logY) 
c 

00 

*  ec  J  e  u  logu  du  +  y. 

0 

Hence 

EXc  -  y  +  eC( Yc  -  y) , 

where  y  -  -  j0  e  u  logu  du,  which  can  be  interpreted  as  the 
C  0 

Euler-Mascheroni  constant  on  [0,c]  (recall  that  y  »  y) .  The  y  can  be 

*  c 

computed  by  numerical  integration;  it  is  a  positive  continuous  function 
in  c  that  increases  on  {0,1]  and  decreases  on  [!,«•]•  Furthermore,  usi 


the  binomial  expansion,  we  have 


E(  exp(  cK  ))  -  E(1  +  c/Y)°-  E  (“)  c®"  T(1  +  k  -  a),  for  a<  1. 
c  k-1  K 


We  end  this  section  with  further  Insights  into  the  irregular 

behavior  of  M  .  Can  M  have  a  limit  distribution  when  p  does  not 
n  n  n 

converge  to  one?  Can  M  have  a  discrete  limit  distribution?  The 

n 

following  results  show  that  the  answer  to  each  of  these  questions  is 


PROPOSITION  3.5.  Suppose  the  birth  and  death  processes  satisfy  (3.1), 

(3.2)  and  p  ♦  p  <  1.  Let  a  be  a  sequence  of  integers.  Then  the 
n  n 


iV»v«»75 


V/W/i w ^ ^ -/a'* .  a • 


distribution  of  -  a^  converges  weakly  to  a  nondegenerate  limit  H  as  n 


•  if  and  only  if 


(3.8) 


In  this  case, 


,  s  -a 

-1  n  n  .  - 

nr  p  +  a  >  0  as  n  +  ». 

ns  n 
n 


H(x)  =  exp(-a  ^(l-p)p^),  xeR, 


where  [x]  denotes  the  integer  part  of  x  (H  is  concentrated  on  the 
integers) . 


Proof.  With  no  loss  in  generality,  we  may  assume  that  p  <  1  for  each  n 

n 


By  (3.6)  with  ®n(x)  =  a^  +  [x],  we  have 


s  -a  -[xj 


n(  1  -  F  (a  +  x))  =  (o(l)  +  n  r  (p 
n  n  ns  n 

n 


-  p)/(l  -  p))  \ 
n  n 


Recall  that  nr  +0  because  of  (3.2).  Then  clearly 
ns 

n 

n(  1  -  F  (a  +x))  ♦  a  *(l-p)p^  if  and  only  if  (3.8)  holds.  Hence,  the 
n  n 

assertions  follow  by  Criterion  2.2. 

Proposition  3.5  may  not  be  too  useful  for  applications  since  (3.8), 

apparently,  is  rarely  satisfied.  Indeed,  we  know  from  Theorem  2.3  ( i) 

that  there  do  not  exist  a  that  satisfy  (3.8)  when  X  ,  ,  p  ,  s  and  p 

n  nk  nk  n  n 

are  Independent  of  n.  A  subtle  variation  in  these  quantities  as  n+«»  is 
therefore  needed  for  the  existence  of  an  that  satisfy  (3.8).  Here  is 
such  an  instance. 

EXAMPLE  3.6.  Consider  the  special  case  in  which  X  ./p  »  n 

nk  nk  * 

s 

k  >  s,  and  II  X  J  p  .  M  >  0  as  n  ♦  «.  Then  (3.1)  and  (3.2)  hold  and 
nl  ni 

Pn*  exp(-logn/[ logn] )  +  e  Let  aR  ■  [logn].  Then 


18 


s-a 

-1  n 

n  r  p 
ns  n 


s  —  i  —  g 

r  p  ♦  Y  e  as  n  + 
ns  n 


Hence  Proposition  3.5  yields 


lim  P(Mn  -  a^<  x)  =  exp(-8e  J),  xeR, 
n 

S  ”1 

where  3  =  ye  (1-e  ).  This  limit  distribution  is  a  discrete  version  of 


G  (x)  »  exp(-e  ). 


4.  Extreme  Values  of  Transient  Processes 


Consider  a  birth  and  death  process,  as  in  Section  2,  with  rates 


A.  ,  u,  .  Assume  that  the  process  is  transient,  that  is,  the  sum  B  *  1  r 

K  k=0  k 

is  finite.  Let  M  denote  the  maximum  of  the  process  up  to  the  time  T  of 
n  n 

its  nth  return  to  state  0.  Because  the  process  is  transient,  the  and 

T  may  be  infinite:  from  (2.1)  we  know  that 
n 

x  -i 

P(M  <  x)  -  I  -  (  £  r  )  ,  and  so  P(M  <  •)  -  1  -  1/B  <  1.  Of  interest, 

k-O 

therefore,  is  the  asymptotic  behavior  of  M  conditioned  on  M  <  °®. 

n  n 

Accordingly,  we  now  consider  the  convergence  of  the  conditional 


distribution 


(4.1) 


where 


P((M  -  a  )/b  <  x  |  M  <  «)  -  H(a  +  b  x)n, 
n  n  n  1  n  n  n 


(4.2)  H(x)  -  P(M  <  x  |  M  <  -)  -  (1  -  (  Z  r  )-1 ) /( 1-1/ B) . 

11  k-o 

Similar  to  the  terminology  above,  we  say  that  conditioned  on 
<  •  has  a  limit  distribution  or  doesn’t  have  one  according  to  whether 
or  not  the  distribution  (4.1)  converges  weakly  to  a  nondegenerate 


distribution.  The  following  result  is  analogous  to  Theorem  2.3.  Here  we 


use  p  -  liminf  V\*  p  =  U^suP  VV  and  8mk  =  11 

k+«  k+® 

THEOREM  4.1.  (i)  The  M  conditioned  on  M  <  «  does  not  have  a  limit 

n  n 

distribution  if  and  only  if 


(4.3) 


liminf  Z  3  ,  <  ». 
,  mk 
m  +<*>  k=m 


CD 

(ii)  Inequality  (4.3)  holds  when  limsup  3^  <  <»  or  when  p  >  1. 

m-K» 

(iii)  Inequality  (4.3)  does  not  hold  when  £.  n  liminf  3  =  •  or  when  p 

k=0  mk 

m-H» 

<  1.  The  condition  p  <  1  is  equivalent  to  lim  *  1. 

k-H» 

(iv)  If  3^  *  lim  3mk  exists  for  each  k,  then  inequality  (4.3)  holds  if 

m  x10 

and  only  if  0^  <  °°. 

Proof.  From  (4.2)  and  a  little  algebra,  we  get 

•  m  « 

(4.4)  (H( x)  -  H(x  -  !))/(!  -  H(x  -  1))  -  (  E  r  /  E  r  )(  E  0  )  . 

k=0  k=0  k=m  m 

Then  arguing  as  in  the  proof  of  Theorem  2.3  (i),  it  follows  that  M 

n 

conditioned  on  M  <  «° does  not  have  a  limit  distribution  if  and  only  if 
n  7 

expression  (4.4)  does  not  converge  to  zero,  which  is  equivalent  to  (4.3) 

(the  first  term  on  the  right  of  (4.4)  converges  to  one). 

Part  (iii)  follows  since,  by  Fatou's  lemma  and  the  product  form  of 

3  .  ,  we  have 
mk 

00  CO  00 

liminf  Z  0  .  >  E  liminf  3  ,  >  Z  b_k, 

ra-*»  k=*m  k*0  m-**  k*0 

and  the  second  assertion  in  (iii)  is  proved  the  same  way  its  analogue  in 


<3 


VMNV 


20 


Theorem  2.3  (iii)  was.  Part  ( ii)  follows  by  a  similar  argument,  and  part 

( iv)  is  a  consequence  of  (ii)  and  (iii). 

Theorem  4.1  says  that  for  a  typical  transient  process  with  liminf 

k+«® 

X,  /  u,  >1,  there  does  not  exist  a  limit  distribution  for  M  conditioned 
k  1  1  n 

on  M  <  ",  but  there  might  be  a  limit  for  a  process  with  lira  X./  =  1. 

n  k+®> 

When  it  is  possible  for  the  limit  to  exist,  then  the  asymptotic  behavior 

of  M  conditioned  on  M  <  »  is  analogous  to  Properties  2.4.  Here  we  have 
n  n 

a  further  simplification. 

REMARK  4.2.  The  distribution  (4.1)  has  the  same  limiting  behavior  as  the 

n  _j  * 

distribution  H(a  +  b  x)  ,  where  H(x)  =  B  E  r.  .  This  follows  since 
n  n  >  n  R 

k=0 

i  x 

1  -  H(x)  =  (1  -  H(x) ) /( ( 1  -  B_1 )  E  r.), 

k=0 

and  so  ratios  of  the  form  (1  -  H(xn>/(1  -  H(yn>),  like  (2.6)  and  (2.7), 
have  the  same  limiting  behavior  as  ( 1  -  H(xn))/(1  -  H(yn)  when  x^ ,  y^  + 


The  preceding  observations  lead  to  the  study  of  when  the  ratio 

Xk/pk  depends  on  n  and  is  nearly  unity  for  large  n.  Accordingly, 

consider  the  maximum  M  ,  as  in  Section  3,  for  a  sequence  of  birth  and 

n 

death  processes  with  rates  X^*  Mn^  that  depend  on  n  as  well  as  k.  Then 


where 


P(M  <  X  |  M  <  •)  «  H  ( x) n 


H  ( x)  -  (1  -  (  E  r  .;)/(!  -  1 / B  ) , 


and  r 


'  %*t  V*  Si"  V 

»  *  (*  »  *  m  w 


j  V  v) j*  o  v  •  .■  ■/  v  v 


"i  '■!"  - :  TV  '-*■ 
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‘  nk 


(u. 


nl 


ynk)/(Xnl  •**  W* 


E  r  .. 
k=0  nk 


Assume  that  (3.1)  and  (3.2)  hold  and  that  p  >  1  for  each  n.  Then  each 

n 

birth  and  death  process  Is  transient  since  B  Is  finite: 

n 


B  =  r  E  p“m  =  r  /(1-1/p  ). 
n  ns  n  ns  n 

n  m=0  n 


(4.5)  B_  =  r  Z 

The  following  result,  analogous  to  Theorems  3.1  and  3.2  combined. 


says  that  when  pn  -*■  1 ,  the  possible  limit  distributions  for 

conditioned  on  M  <  »  are  G  ,  0  <  c  <  ».  Here  we  let 
n  c 


c  =n(l-l/p)r* 
n  n  ns 


c  =  1  im  c 


n 


and 


a  =  s  -  1 
n  n 


a  =  s  -  1 
n  n 


an  =  sn  -  1  -  log(l/cn  -  l/n)/logp 


,  -1 
b  =  nr 
n  ns 


bn  =  l/logpn 


bn  =  l/logpn 


when  c  =  0 

when  0  <  c  <  <» 

when  c  =  °°. 


E3 


yy. 

L" 

m 


•  * . 

y; ' 

1 

& 

fcf; 
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This  notation  is  that  of  Section  3  with  p^  replaced  by  1 /p^  and  the  last 

a  changed  slightly, 
n 

THEOREM  4.3.  Suppose  (3.1)  and  (3.2)  hold,  p^  >  1  for  each  n,  and 

P  +1.  Then  M  conditioned  on  M  <  »  has  a  limit  distribution  if  and 
n  n  n 

only  if  c  +  c  where  0  <  c  <  «.  In  this  case 
n 


(4.6) 


lim  P(M  -  a  )/b  <  x  I  M  <  “)  =  G  (x)  ,  xeR, 

n  n  n  n  c 

n  -k® 


where  a  ,b  are  defined  above, 
n’  n 

Proof .  We  will  prove  that  c^  ♦  c  implies  (4.6).  Then  the  rest  of  the 
assertion  will  follow  by  the  argument  we  used  in  the  proof  of  Theorem 


n 


3.2. 

Similar  to 

(4.5),  we  can  write, 

m 

m-s 

n  » 

(4.7) 

E  r  . 

1  n  nk 

k*0 

=  r  Z  p  ■  r 

nsn*=0  n 

-m+s  -1 

(1  -  Pn  °  )/( 1  -  l/pn). 

n 


Let  ran(x)  denote  the  integer  part  of  +  b^x.  Then  using  (4.5)  and 


(4.7),  we  hav  e 


m  (x) 
n 


(4.8)  n(l  -  Hn(an  +  bRx))  =  n(Bn/  Z  rnk  -  1)/(B  -1) 

k=0 

-m  (x)+s  -1 

=  Cn((1  “  pn  n  n  }  '1)/(1  "  cn/n) 

m  (x)-s  +1 

=  l/((l/c  -  1/n) ( p  n  n  -  1)). 

n  n 

Suppose  c  c,  where  0  <  c  <  «.  Then  by  Cases  1  and  2  in  the  proof 

of  Theorem  3.1,  with  p  replaced  by  p  *,  it  follows  that 

n  n 

m  (x)-s 

(4.9)  lim  n(l  -  H  (a  +  b  x))  *  lim  c  /(p  n  n-l) 

n  n  n  n 

n  +<*•  n 

=»  x  ^  when  c  =  0 

X 

=  c/(e  -  1)  when  0  <  c  <  ®. 

Next,  suppose  +  «°.  Here 

m  (x)  =  s  -  1  -  log(  1  / c  -  l/n)/logp  +  x/logp  +  0(1). 
n  n  n  n  n 


Using  this  in  (4.8),  we  have 


m  (x)-s  +1 
n  n 


(4.10)  lim  n(l  -  H  (a  +  b  x))  =  lira  l/(l/c  -l/n)p 
,  n  n  n  n  n 

n+“ 

-  lim  exp(-(m  (x)  -  s  +  l)logp  -  log(l/c  -  1/n)) 
n  n  n  n 

n-H* 


Thus,  (4.9),  (4.10)  and  Criterion  2.2  yield  (4.6). 

REMARK  4.4.  The  analogue  of  Approximation  3.3  for  the  birth  and  death 

process  with  X^/ =  p,  k  >  s,  and  p  >  1  is  as  follows 

(4.11)  P((Mn  -  s  +  l)logpn  <  x  |  Mn  <  »)  s  exp(cn/(ex  -  1)),  x  >0, 


where  c  =  n(l  -  l/p)  n  A./u  , 

k=l  *  * 


5.  Extreme  Values  of  Queues 


We  now  apply  the  preceding  results  to  the  M/M/s  and  related  queues. 
The  M/M/s  queueing  process  described  in  Section  1  is  a  birth  and  death 
process  with  birth  rate  X^  =  X  (the  Poisson  arrival  rate  of  customers) 
and  death  rate  p^  =  p  min{k,s}  (the  rate  at  which  k  customers  depart  from 
the  s  servers). 

For  our  first  result,  we  suppose  M  is  the  maximum  of  this  M/M/s 

n 

queue  up  to  the  nth  time  the  system  becomes  empty.  The  limiting  behavior 
of  M^  depends  on  the  queue's  traffic  intensity  p  =  X/sp.  The  queueing 
process  is  positive  recurrent  when  p  <  1,  null  recurrent  when  p  =  1,  and 
transient  when  p  >  1.  The  following  is  an  immediate  consequence  of 
Theorems  2.3,  2.5  and  4.2. 

COROLLARY  5.1.  If  p  <  1 ,  then  M  does  not  have  a  limit  distribution.  If 

n 

p  =  1 ,  then 


(5.1) 


_  -1 

lira  P(mn/nb  <  x)  =  e  X  ,  x  >  0, 
n  +«>  1 


g 

where  b  =  (X/p)  s!.  If  p  >  1 ,  then  M  conditioned  on  M  <  ®  does  not 

n  n 

have  a  limit  distribution. 

For  the  next  result,  we  suppose  that  Mn,  as  in  Sections  3  and  4,  is 

the  maximum  for  an  M/M/s  queue  with  arrival  rate  X  ,  =  X( n) ,  service 

n  nk 

rate  p  .  =  p(n)min{k,s  },  and  number  of  servers  s  .  The  traffic 
nk  n  n 

intensity  of  the  queue  is  p  *  X(n)/(s  p(n)).  Clearly  X  ,/p  .  =  p  .  for 

n  n  nk  nk  n 

k  >  s  .  We  will  use 
n 


COROLLARY  5.2. 


ns 


if 

*■1 


Suppose  n 


-1 


uni/XrU 


s 

n 


r  +0  and 
ns 

n 


!  (  p(n)/  X(n)  ) 

p  *  1 . 
n 


s 

n 


(i)  If  p  <  1  for  each  n,  then  M  has  a  limit  distribution  if  and  only 

n  ’  n  J 

if  n(l  -  p  )r  *  -*•  c,  where  0  <  c  <  °°.  In  this  case,  the  limit  distri- 

n  ns 

n 

bution  is  and  appropriate  norming  constants  are  as  in  Theorem  3.1. 

(ii)  If  p  >  1  for  each  n,  then  M  conditioned  on  M  <  <*>  has  a  limit 

n  n  n 

distribution  if  and  only  if  n(  1  -  1/p  )r  *  c,  where  0  <  c  <  «>.  In 

n  ns 

n 

this  case,  the  limit  distribution  is  G and  appropriate  norming  constants 
are  as  in  Theorm  4.3. 

Proof .  The  two  assertions  are  special  cases  of  Theorems  3.1,  3.2  and 
Theorem  4.3,  respectively.  Note  that  condition  (3.1)  is  satisfied,  and 
so  is  (3.2)  since 

sn 

n  *  E  r  .  <  2n  *r  -*•  0. 

,  n  nk  ns 

k=0  n 

REMARKS  5.3.  (a)  The  number  of  customers  in  an  M/M/°°  queueing  system 

over  time  is  a  birth  and  death  process  with  rates  X^  *  X  (the  Poisson 

arrival  rate  of  customers)  and  p^  =  kp  (where  p  is  the  service  rate  of 

each  of  the  infinite  servers).  The  traffic  intensity  is  p  =  X/ p.  The 

first  and  third  assertions  of  Corollary  5.1  also  hold  for  this  queue,  but 

there  are  apparently  no  analogues  of  (5.1)  or  Corollary  5.2. 

b)  Consider  a  service  system  in  which  the  number  of  customers  in  the 

system  over  time  is  a  birth  and  death  process  with  rates  X^,  p^  that 

represent  customer  arrival  and  departure  rates  when  k  customers  are 

present.  We  assume  that  X  /  p.  =  P  for  k  >  s,  where  s  is  a  specific 

p  k 

state,  but  we  place  no  other  restriction  on  the  rates.  We  call  this  an 


M/M/GR-s  queueing  process,  where  GR  stands  for  general  rates.  General 
rates  are  used  for  modeling  such  phenomena  as  balking  and  reneging  of 
customers;  non-standard  service  disciplines;  dynamically  changing  rates 


under  a  control  policy  that  minimizes  the  system  cost;  and  simultaneous 

customer  processing,  where  is  the  total  workrate  when  1c  customers  are 

present  and  customer  i  receives  p  of  the  workrate  (p,+...+p  =  1). 

1  1  m 

Corollaries  5.1  and  5.2  readily  extend  to  M/M/GR-s  queues. 

(c)  The  approximations  (3.7)  and  (4.11)  apply  to  the  M/M/s  and  M/M/GR-s 
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